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ABSTRACT 


A study of two oceanic extratropical cyclones, generated 
by the Navy Operational tmospheric Prediction System 
(NOGAPS) model (see Sandgathe, 1981), ravaals a simultaneous 
development of the upper- and lower-tropospheric ‘features 
while these features maintain their positions relative to 
each other. These cyclones also develop extreme shear in 
most synoptic fields after 120 hours of numerical 
Simulation. 

The mass budgets of these cyclones show a concentration 
of inward lateral transport (convergence) in the lowest 
model layer and outward lateral transport (divergence) ina 
layer centered near 300 mb. Tine sections of lateral trans- 
port, vertical velocity and mass tendency reveal that these 
cyclones develop in two phases--an explosive cyclogenesis 
phase, and a quasi-steady state phase. These time sections 
also reveal a 12-hour cyclic pattern embeided within the two 
phases of cyclone evolution. This 12-hour cyclic pattern is 


not evident in the map sequences of the surface pressure. 
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I. INTRODUCTION 


For many years maritime cyclones have been of interest 
to military and civilian communities from both a scientific 
and a practical point of view. Oceanic lows frequently 
develop rapidly and produce meteorological phenomena which 
can adversely impact civilian and military operations in 
oceanic and coastal regions. Accurate and dependable pre- 
diction of the development and movement of oceanic cyclones 
is required for reliable maritine forecasts. 

Numerous studies of middle-latitude cyclones have 
focused on continental cyclones (e.g. Smith (1973), 
Petterssen and Smebye (1971), Johnson and Downey (1976), 
Chen and Bosart (1977), and others). This 1S primarily due 
to the abundant observational jiata which is readily avail- 
able from continental regions. In contrast, relatively few 
studies have investigated maritime cyclones because of the 
limited observational data from oceanic regions. As a 
result, the factors influencing maritimS cyclogenesis are 
not fully understood. 

Several studies of oceanic cyclones indicate that sig- 
nificant differences exist in the physics involved in conti- 


nental and oceanic cyclogenesis. Petterssen, Bradbury and 
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Pedersen (1962) compared the development of cyclones over 
the western North Atlantic to that of cyclones over the 
eastern United States. This observational study indicates 
that the initial development of oceanic cyclones appeared +o 
result primarily from low-level thermal advection while con- 
tinental cyclones usually developed as a result of an 
upper-level vorticity maximum approaching a low-level baro- 
clinic zone. A study by Nitta and Yamamoto (1974) also sug- 
gests that maritime cyclones frequently develop in the 
vicinity of an ocean baroclinic zone in conjunction with 
nearly zonal upper-level flow. Sanders and Gyakum (1980) 
found that explosive cyclogenesis 1s preiominantly a mari- 
time, cold-season event and suggested that the roles played 
by cumulus convection and air-sea fluxes in cyclogenesis are 
not adequately understood. A numerical study by Sandgathe 
(1981) indicates that air-sea fluxes play a Significant role 
in the evolution of oceanic cyclones. These studies and 
others suggest that the air-sea fluxes of latent and sensi- 
ble — play a Significant role in oceanic storm 


development. 
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Current operational numerical models have been shown to 
be seriously deficient in predicting maritime cyclone 
deepening rates. The Sanders and Gyakum (1980) study indi-~ 
cated that the National Meteorological Canter (NMC) primi- 
tive-equation models predicted less than one half of the 
deepening rates observed in explosive oc2zanic cyclogenesis. 
The failure of numerical models to predict explosive cyclo-~ 
genesis is not entirely the result of the limited observa- 
tional data from oceanic regions. Part of the problem lies 
in modeling vertical exchanges in the atmosphere. The upper 
troposphere is dominated by wave structures, while vortex 
structures dominate in the lower troposphere. Most cases of 
upper tropospheric wave amplification are accompanied by the 
development of anticyclones and cyclones in the lower tropo- 
sphere. This suggests that th2re exists a strong interac~ 
tion between these two regimes. The smaller size of oceanic 
cyclones presents another probl2m. The space scale of upper 
tropospheric waves is typically two to four times the space 
scale of oceanic cyclones. As a result, it appears chat 
numerical models adequately Simulate the large-scale 
processes involved in synoptic-scale systems, but do not 
capture the smaller-scale features embedded within 


synoptic-scale weather systems. 
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Due to the sparse observational network in oceanic 
regicns, few quantitative studies of maritime cyclogenesis 
have been completed. To fill this void, studies of numeri- 
cally simulated cyclone development are needed. Numerical 
models, having an internally consistent "Jata" base, elimi- 
nate errors resulting from incomplete or inaccurate observa- 
tions from data-sparse regions such as dceans. In cases 
where there is a favorable comparison between numerically 
simulated cyclogenesis and real atmospher2 cyclogenesis, the 
model may provide a greater understanding of the physics 
involved in cyclone development than would the incomplete 
observations of the real atmosphere storm. A disagreement 
between the model and the real atmospher= would not neces- 
Sarily prove to be of little or no value. In cases where 
the model and the real atmosphere do not agree, the diaanos- 
tics of the real atmosphere may well isolate the reasons for 
the disagreement and thus lead to improvement in the model. 

A basic feature of cyclogenesis iS an inward mass trans- 
port (convergence) in the lower troposphere and and outward 
mass transport (divergence) in the upper troposphere. Asso- 
ciated with these horizontal mass transports is an upward 


vertical mass transport within ‘the cyclone. Quantifying 
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this mass circulation can provide important information about 
the large-scale forcing within 2 cyclone. 

Cyclones are vortices embedded in a planetary-scale 
pattern. Consequently, a study of the basic properties of 
cyclenes should take advantage of their vortex nature. One 
such approach is the investigation of the budgets of the 
fundamental quantities such as mass and circulation within a 
volume moving with the storm. Quasi-Lagrangian transport 
theory, outlined py Johnson and Downey (1975 a and b), is 
particularly suited for the stuiy of cyclones. 

In this thesis the mass budgets of nunerically simulated 
eceanic cyclones are investigated to ascertain the basic 
aspects of cyclone development and maintenance. A quanti- 
tative study of the mass circulation helps describe the 
ihtensity and physical structure of the cyclone and Weta ce 
vertical variations in the cyclone circulation. 

This study follows the Ph.D. thesis investigation by 
Sandgathe (1981). The budget studies were applied to two 
storms generated by a version of the UCLA general circula- 
tion model, which includes a sophisticated parameterization 
scheme for planetary boundary layer (PBL) and free atmo- 


sphere interactions. The Sandgathe study demonstrated the 
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role of air-sea fluxes of heat and moisture in initial 
cyclogenesis. He showed that cyclon2 intensification 
resulted when sensible heat fluxes were removed. In con- 
east, the lack of moisture fluxes rasulted in cyclone 
decay. 

In this thesis, two cyclones were studied--one in the 
southern (spring) hemisphere and one in the northern (fall) 
hemisphere. Both storms are characterized by rapid 
development during the period from 96 to 126 hours of simu- 
lated cyclogenesis. In each case, the central surface pres- 
sure of these storms decreases at a rat2 of nearly one mb 
per hour during this 30 hour period. Thus, these storms are 
good examples of the explosive deepening which is not 
adequately handled by current operational dynamic models. 

This study is one in a series of experiments designed to 
investigate *he physical processes involved in extratropical 
cyclone development. The primary objectives of this study 
are: 


e +¢+o study the synoptic features of two model-generated 
cyclones and compare them to observed cyclones; 


e +o isolate the lateral and vertical branches of the 
Mass circulation elibherat jelay different stages in the 
evolution of a model-generated maritime cyclone; and, 


e +¢o compare vertical velocity values derived through 
quasi-Lagrangian diagnostics with the vertical velocit 
values rom the modéls internally consistent "data 
set. 
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The atmospheric model used for this study is briefly 
described in Chapter II. The results of this study begin in 
Chapter III with an objective review of the synoptic fea- 
tures of the two model storms. A description of the quasi- 
Lagrangian diagnostics technique and the equations used in 
the mass and circulation budgets are presented in Chapter 
Iv. The generalized budget equation and the budget volume 
coordinate system are described in Appendix B. Chapter V 
contains the analysis of the mass budget results. Cone ius 
sions and recommendations for further studies are presented 


in Chapters VI and VII respectively. 
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Il. ATMOSPHERIC 42DEL DESCRIPTION 


$< MODEL STRUCTURE 

The Sandgathe experiments used the Navy Operational 
Atmospheric Prediction System (NOGAPS) model which is being 
developed by the Naval Environmental Prediction Research 
Facility (NEPRF). The NOGAPS model is a state-of-the-art 
version of the UCLA general circulation model. A brief 
overview of the NOGAPS model is presentei in this chapter, 
Por aaore detailed description of the NOGAPS model the 
reader is directed to Sandgath2 (1981). The UCLA model is 
d2scribed in detail by Arakawa and Lamb (1977). 

A fine resolution grid (2 degrees latitude by 2.5 
degrees longitude) within a 60 degree longitude global 
sector is used in this model. This grid resolution, which 
equates to a 222 km by 196 km grid distance at 45 degrees 
latitude, is sufficient to resolve shorter cyclone waves. 

The prognostic variables of the molel include the sur- 
face pr2ssure minus 50 ob, 7, the horizontal wind velocity, 
V, temperature, T, and specific humidity, q. Several other 
prognostic variables (to be discussed lat2r) are associated 


with the PBL. 


Ws, 





The model uses a sigma vertical coordinate system 


defined as: 





P= P 
= 2 

T 

where: 
P = pressure 
ie 
P. = surface pressure 
P. = pressure at the top of the model atmosphere (50 mb). 


Thus the top and bottomy moundaries of the model area o = 0 
and 9 = 1. The NOGAPS model consists of six vertical layers 
(Fig. 1). In this depiction, an arbitrary surface pressure 
of 1000 mb has been assumed. Variables carried at the layer 
mid-levwels are T, V, and q. The vertical velocity, io 2 Aue 
defined at the layer interfaces. 

The model employes a staggered grid (Arakawa scheme C) 
im” the horizontal (Fig. 2). The variables defined at the 
Genter polints are T, gq, and WO. The meridional wind compo- 
nent, v, is carried one degree north and south of the center 
point with «he zonal wind component, u, carried 1.25 degrees 
Sesteapa west of the center point. For diabatic and fric- 
tional computations, the u and v components are averaged to 
+he center pcint. A smoothing of the meridional component 


adjacent to the poles is implemented in theis nodel. 


— 
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The area used for this study, 60 degrees wide extending 
from pole to pole, is an open-ocean (i.e. entirely over 
Water with no sea ice) domain which is fixed in time with 
two zones (centered at 45N and 35S) of maximum north-south 
gradients in sea-surface temperature. Continuous over-ocean 
propagation within the model is guaranteed by imposing 
cyclic boundary conditions in the east-west direction. 
There are 91 and 24 grid points in the north-south and 
east-west directions respectively. For this thesis, zero 
degrees longitude has arbitrarily been assigned to the 
western boundary. 

The PBL paramererization scheme of the NOGAPS model, 
Which follows Deardorff& (1972) and Randall (1976), provides 
for interactions between the PBL and cumulus and/or stratus 
cloud ensembles at each grid point. Momentum, heat and 
moisture fluxes at the surface are diagnostically determined 
uSing the values of T, ¥V, andq from the dynamic portion of 
*he model anda bulk Richardson number which is based no 
sea-surface temperature. The cumulus parameterization 
scheme follows Arakawa and Schubert (1974) and Lord (1978). 
Me®ecees Of water vapor, predicted cloud distributions, and 
large-scale precipitation processes are included in the 


radiational heating package which follows Katayama (1972). 
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Figure 2. Horizontal distribution of model large-scale 
prognostic variables. This distribution follows 
Arakawa scheme C. (After Sandgathe, 1981) 


B. INITIAL CONDITIONS 

The Sandgathe study utilized initial conditions which 
are representative of open-ocean fall and spring regimes. 
foe Gitta! conditions for the model atmosphere follow 


Simmons and Hoskins (1977, 1978). Wind PEOrcILecs Fice 


the fall hemisphere were determined using a polynomial 
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approximation of the January mean zonal wind profile at 30N 
from Oort and Rasmussen (1971). The fall hemisphere winds 
are multiplied by a factor of 1.5 at each level to provide 
the spring hemisphere wind profiles. 

The model start time was 1 IJctober. Thus, with respect 
to solar fluxes, the northern and southern hemispheres are 
the fall and spring regimes respectively. Sea~surface 
temperatures and atmospheric jets are also representative of 
these seasons. Zonal jets of 35 om secu~!? and 52 mm secn!, 
centered at 45N and 35S respectively, represent possible 
conditions for these regimes. As shown by Simmons and 
Hoskins (1977, 1978), the existance of these jets guarantees 
the growth of baroclinically unstable waves. 

Sandgathe used the conditions listei above as initial 
conditions in his experiments on the evolution of extratrop- 
ical cyclones. A small amplitude sinusoidal perturbation in 
wave number six was added to the meridional wind component 
at each level. The maximum value of the perturbation was .5 
msecm! at 45N/S. The perturbation at the equator and the 
poles was zero. The associated perturbation in the surface 


pressure field was approximately one millibar. 
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III. SYNOPTIC DISCUSSION 


A. OVER VIEW 

An initial objective of this thesis is a comparison of 

the synoptic features of numerically simulated and observed 

ceanic cyclones. This discussion will focus on the signif- 
icant features in the model "data" fields. Available "data" 
consists of model output at three-hour intervals from hour 
96 to hour 168 inclusive. Except for the 1000 mb heights 
and surface pressures, all fields in this discussion are 
generated from mid-layer model data. With the exception of 
+he wind fields, a three point zonal smoothing technique was 
applied to all data fields represented in this chapter. The 
u and v wind components were averaged to the height grid 
points for the wind field representations. 

The synoptic summary of th2 two model cyclones begins 
With a presentation of the significant features in the lower 
troposphere using surface pressure and temperature fields, 
the 1000 mb height fields and the 1000 mb to mid-layer 4 
(approximately 1000-500 mb) thickness fields (Figs. 4 to 


11). This is followed by a summary of mid-tropospheric 
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features (mid-layer 4, approximately 590 mb) of the height, 
temperature, wind, and vorticity fields (Figs. 12 to 19). 
Significant features in the mid-layer 3 (approximately 300 
mb) height, wind, and temperature fields are described in 
the fourth section (Pigs. 20 to 27). This chapter concludes 
with a brief comparison of the development of the model 
cycicnes and observed cyclones. 

Table I contains a list of synoptis fields relevent to 
this summary. Fields valid at 96, 120, 144, and 168 hours 
after the initial time are presented. With respect to the 
Synoptic summary and the figures presented in this thesis, 
the indicated times represent the number of hours of numer- 
ical simulation. The storm tracks of the two numerically 
‘Simulated cyclones used for the budget studies of this 
thesis are shown in Fig. 3. 

The figures presented in this chapter depict two 


adjacent 60 degree sectors bas2d on the cyclic continuity 


assumption of the model. Six and nine degree budget volume 
Circles, centered on the surface pr2ssure mininuna, are 
Gepicted om each chart. Radii six and nine were chosen 


because they include the more intenss portion of the 


cyclenes during the early and latter stages cf the study 
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respectively. In following discussions, inner radius and 
outer radius refer to the six and nine degree budget radii 
respectively. Both the northern hemisphere (NH) and south- 
ern hemisphere (SH) cyclones are discussed, relative +o each 
figure, to point out differences and similarities in their 


development. 


TABLE [I 
DATA FIELDS DEPICTED IN THE SYNOPTIC SUMMARY 


SIGMA APPROXIMATE 
LEVEL PRESSURE 
(OR LAYER) LEVEL(S) FIELDS 
level 6 surface pressure (mb) / 
temperature 
(1000 mb)/ 1000/ 1000 - 500 mb 
mid-layer 4 500 mab thickness/ . 
1000 mb height 
level 6 surface pressure/ 
temperature 
mid-layer 4 500 mb WOR CA Car y 7, 
height 
mid-layer 4 500 mb temperature/ 
height 
1d- > 300 mb temperature 
mid-layer 3 m eavcne / 
mid-layer 3 300 mb wind vectors/ 
isotachs 
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Figure 3. Storm tracks of two model cyclones. Indicated 


a= 8 


mid- 


locations are for six hour intervals from hour 96 
Poecowgn nour 168. Crrcle radius is 6° latitude. 
a. Northern hemisphere. b. Southern hemisphere. 


LOWER—=TROPOSPHERE FEATURES 
Surface pressure and temperature fields valid at hour 96 
sheovm in Fig. 4. The 1000 mb height and 1000 mb to 


layer 4 (approximately 1000 - 500 mb) *hickness fields 
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for the same time are presented in Fig. 5. Weak thermal 
ridges, at 32 (92) and 57 (117) degrees longitude respec- 
tively, are apparent in the NH. A 1011 mb low (cyclone N1) 
is lecated in the thermal ridge near 32 (92) degrees longi- 
tude. A somewhat more intense Low (cyclone N2), with a cen- 
tral surface pressure of 1007 mb, is in the thermal ridge 
near 57 (117) degrees longitude. rnitially, cyclone N2 
appears to be the main feature; however, within 24 hours, 
cyclone N1 becomes the dominant low in the NH sector. A 
third area of low pressure approximately 1000 nm northeast 
of cyclone N1 develops very little, but does persist through 
hour 156. This Low pressure area will be discussed briefly 
later in this section. In the SH, a 1003 nb low (cyclone S) 
is located in a thermal ridge oriented south-southwest (SSW) 
to north-northeast (NNE) and centered near 40 degrees longi- 
tude. Table II contains the central surface pressure values 
of the three cyclones discussei in this synoptic summary. 
Cyclcnes hl, N2, and S$ all display typical lower- 
tropospheric thermal advection patterns with warm advection 


ahead of, and cold advection behind and equatorwardof the 


developing cyclones. 
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TABLE If 


CENTRAL SURFACE PRESSURES OF THREE MODEL CYCLONES 


(in mb) 
= CYCLONE CYCLONE CYCLONE 
96 1011 1007 1003 
i 1008 1002 1001 
10 2 1006 1001 998 
105 1006 1002 996 
108 1005 1001 994 
111 1000 998 oe 
114 997 996 280 
117 995 998 es 
128 932 998 986 
W2.3 988 997 984 
126 987 997 or 
w2-9 988 998 oS 
132 989 1000 O79 
135 989 1001 978 
13 8 987 1001 976 
14 1 ao 1004 975 
144 a0 1005 975 
14 7 988 1006 975 
150 987 1006 976 
15 3 988 1008 978 
156 988 1008 980 
ipo 985 1006 978 
16 2 98/5 1006 o) Is) 
165 985 1006 980 
16 8 984 1006 981 
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Surface pressure and temperature fields at hour 120 are 
given in Fig. 6 with height and thickness fields in Fig. 7. 
During the period from hour 96 to ho ur 12), the NH lows move 
appreximately 700 nm to the east while moving only slightly 
poleward. Cycione N1 decreases to 992 ab and becomes the 
Main surface low, while cyclone N2 decreases to only 998 mb. 
Approximately 70% of the surface pressure decrease experi- 
enced by cyclone N1 occurs during the latter half of this 24 
hour period. In contrast, more than half of the surface 
pressure decrease experienced by cyclone N2 occurs during 
*he first three hours of this period. It should be noted 
that cyclone N1 experiences 32 rapid central surface pressure 
decrease of 17 mb between hour 108 and hour 123. In the SH, 
cyclone S decreases to 986 mb at a nearly constant rate of 
appreximately two millibars per three hours. 

The thermal ridges associated with cyclones N1 and S are 
more pronounced than at hour 96; however, thermal advection 
at the surface increases very little during «he period from 
nour 96 to hour 120. By the end of this period, pockets of 
warm air are coincident with the low pressure centers in the 
NH. However, cyclone N2 does’ not continue to develop and 


thus will not be discussed further. 
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The analyses for hour 144 are presented in Pigs. 8 and 
9. During the period from hour 120 to hour 144, cyclone N1 
moves approximately 700 nm to the east-northeast (ENE) while 
cyclcne S$ moves approximately 500 nm to the east-southeast 
(ESE). Cyclone N1 experiences little intensification as the 
central surface fluctuates between 992 and 987 mb with the 
lowest value cccurring at hours 126 and 138. The low pres- 
sure area, initially to the northeast of cyclone N1, inten- 
sifies slightly and moves to a location approximatey 1000 no 
to the northwest of cyclone N1 in an area of lower-tropos- 
pheric warm advection. During this same period, the central 
surface pressure of cyclone S decreases to 975 mb with more 
than half of this decrease ocsurring during the first three 
hours of the period. At hour 132 (not shown), a pocket of 
warm air is coincident with the surface center of cyclone S. 
A significant feature of the NH thickness field is the north 
to south orientation and the amplitude of the thermal ridge. 
In contrast, the SH thermal ridge is oriented nearly west to 
east and is extremely elongated. 

Surface pressur2 and *empérature fields valid at hour 
168 are presented “in Fig. 10 with the height and thickness 


fields given in Fig. 11. During the period from hour i144 to 
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hour 156 hours, cyclone N1 remains relatively steady-state 
with a central surface pressure of approximately 988 amb 
while the central surface pressure of cyslone S$ rises 5 mb 
to 980 mb. From hour 156 to hour 168, the central surface 
pressure of cyclone N1 decreases to 984 ab while cyclone §S 
remains nearly constant at about 980 mb. During the final 
24 heurs of this study, the NH thermal ridge shifts from a 
north-south orientation to a northwest-southeast orienta- 
tion. However, the SH thermal ridge maintains its west-east 
orientation. The 1000 mb height contours nearly parallel 
the thickness contours poleward of each cyclone so there is 
little warm advection in the vicinity of these cyclones; 
however, cold air advection continues equatorward of each 
cyclone. During the period from hour 144 to hour 162, 
cyclone N1 moves approximately 300 nm to the NNE. ih 
appears that cyclone N1 moves approximately 360 am during 
the last six hours of this study; however, this apparent 
rapid movement may actually ba a reformation of the low in 
“he vicinity of «he low pressure area Il 5dcated northwest of 
cyclene N1 at hour 144. Cyclone S also exhibits an apparent 
rapid movement during the final 24 hours of simulation, mov- 


ing approximately 720 nm to the ESE. Approximately 500 no 


Sys, 





of this movement occurs during the perioi from hour 150 to 
hour 162 as the central surface pressure increases three 
millibars to 979 mb. As in the NH, it is possible that this 
apparent rapid movement is a reformation sf cyclone S. How- 
ever, unlike the NH case, there is no secondary low or warm 


pocket in the vicinity of the naw location of cyclone S. 


on MID- LEVEL FEATURES 

In this section mid-tropospheric features are discussed 
using mid-layer 4 of the model (approximately 500 mb). 
Mid-layer & height and temperature fields valid at hour 96 
are depicted in Fig. 12 with the vorticity and wind fields 
Given ia Piq. 13. A weak short wave within nearly zonal 
flow is apparent in each ess: tere, The height and t¢emper- 
ature fields are nearly in phase so there is little thermal 
advection in either hemispheres. Developing cyclones N1 and 
S are located under the short wave ridges in the respective 
hemispheres with a weak vorticity maximum located 300 to 400 
mm to the west of each cyclone. Mid-tropospheric wind flow 
patterns indicate weak cyclonic vorticity advection (CVA) 


over each cyclogenesis area. 
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Mid-layer 4 fields valid at hour 120 are presented in 
Pegs. WW and 15. The height and thernaal fields remain 
nearly in phase with weak cold advection immediately to the 
west of the center of cyclone N1 and warm advection approxi- 
mately 5 degrees longitude further to the west. During the 
first 24 hours of this study, the amplitude of the NH short 
wave doubles while that of the SH short wave increases by 
approximately 70% This amplitude intensification is nearly 
Simultaneous with *+he most rapid surface development which 
occurs during the period from hour 96 to hour 126. 

The orientation of the NH Short wave trough is northwest 
+o southeast while that of the SH short wave is southwest to 
northeast as each cyclone Maintains its position relative 
its associated short wave. The cyclonic vorticity maxinun 
in the NH increases in magnitude andis located directly 
over cyclone N1 resulting in anticyclonic vorticity advec- 
tion (AVA) west of the center of cyclone Ni. The increase 
in the vorticity maximum at mid-layer 4 and the most rapid 
surface development also occur Simultaneously. In the SH, 
the cyclonic vorticity maximum also increases in magnitude; 
however, in contrast to the NH case, the SH vorticity max- 


imum does not move over the storm center. Tnstead the 
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vorticity maximum moves to the east as cyclone $ moves 
southeasterly to a position poleward of the vorticity max- 
imum. As in the NH case, intensification of the vorticity 
maximum occurs simultaneously with the most rapid surface 
development. Zero vorticity advection occurs over cyclone §S 
while strong CVA occurs to the north and northeast. 
Mid-tropospheric fields valid at hour 144 are given in 
Figs. 16 and 17. The amplitude of each short wave continues 
to increase and the orientation of the trough axes becomes 
more west-east. Extremely sheared flow 1s evident in the SH 
and, to a lesser deqreé, in the NH. The short waves and 
their associated cyclones maintain their positions relative 
*9 each other during the period between hours 120 and 144 as 
os 
the height and thermal fields remain nearly in phase. 
Although maintaining itS position over cyclone N1, the vor- 
icity maximum in the NH weakens, having reached a peak at 
hour 132 (not shown). As in the NH case, the SH vorticity 
maximum reaches its peak at hour 132 and naintains its posi- 
+ion relative to cyclone 5S. In each case, the vorticity 


patterns are quite elongated, because of the extreme shear 


present in the wind field. 
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At hour 168 (Pigs. 18 and 19), cyclone N1 is positioned 
under the inflection point ahead of the trough. The NH 
short wave continues to amplify and th2 shearing is more 
pronounced as the axis of the short wave is nearly east +o 
west. Between hours 144 and 168, the SH upper trough fills 
leaving a split flow pattern with nearly zonal flow to the 
south and the weakening trough to the north of 46S, the 
det itude of cyclone S. A weak closed low appears approxi- 
mately 700 nm to the northwest, however cyclone S fails to 
develop into a deep cold core vortex. The vorticity fields 
continue to weaken and become more elongated during the 
final 28 hours of this study. Broad, relatively uniforn 
areas of cyclonic vorticity ar2 located squatorward of each 
cyclone. These vorticity patterns nearly parallel the wind 
fiow iadzecating little vorticity advection in either 


hemisphere. 


De UPPER-LEVEL FEATURES 

Mid-layer 3 (approximately 300 mb) height and tempera- 
ture fietd= valid at hour 96 are given in Fig. 20 with the 
wand fields given in Fig. 21. Nearly zonal flow is evident 
with weak cyclonic flow 700 nm to the west of each cyclone 


and weak anticyclonic flow 1000 nm and 2900 nm to the east 
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of cyclones Ni and S respectively. The two jets, centered 
near 45N and 35S, are nearly zonal and somewhat stronger 
than those in the initial field. The NH and SH jet maxima, 
approximately 45 m sec! and 55 m secu! respectively, are 
located 900 to 1500 nm to the west sf the developing 
cyclones. 

During the first 24 hours of this study, the flow imme- 
diately to the west of each cyclone becomes increasingly 
cyclonic as the amplitude 9f the trough/ridge systen 
increases. At hour 120 (Figs. 22 and 23) the cyclones are 
positioned slightly east and poleward of «the base of the 
trough. The jet maximum in the NH increases only slightly 
while the SH jet maximum increases to over 60m _ sec~!. 
These jets are located approximately 6090 om west of the base 
of the respective troughs. 

The upper-level charts valid at 144 hours are shown in 
Pigs. 24 and 25. In the NH, cyclone N1 is located ahead of 
the trough under nearly meridional flow with a jet maximun 
of cver 40 om sec~! located approximately 500 am to the 
southeast and 2 secondary jet maximum, also greater than 40 
msecmt, located near the top of the ridge, approximately 


600 me to th® north of cyclone N1. In the SH, there are two 
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jet maxima (each greater than 50 m sec!) associated with 
cyclone S--one behind the trough and one in the base of the 
trough. In each hemisphere the flow ahead of the trough is 
basically meridional and the short waves are continuing to 
increase in amplitude. 

The mid-layer 3 wind field valid at hour 156 (not shown) 
is characterized by several small disjointed jet maxima 
greater than 40 m secm! in the NH. There appears to be no 
organization to these areas as some are ahead and some 
behind the trough. By this time, the SH trough has tilted 
back on itself and is oriented nearly west to east and the 
flow south of the trough is now zonal. A developing cut-orf 
low is indicated by the extrame shear in the SH height 
field. Similar directional shear is present in the NH, alt- 
hough not to the extent of that in the SH. 

Figs. 26 and 27 depict the mid-layer 3 patterns at hour 
168. In the NH, the position 29f cyclone N1 relative to the 
+rough remains unchanged during the final 24 hours of tais 
study. However, an organized jet maximum of over 40 a secu! 
to the west is approaching th= base of the trough. This 
maximum is a result of the cyclic continuity assumption 


emplcoyed by the model and is the same maximum which was in 
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the downstream ridge (north of the cyclons) at approximately 
hour 144, The flow directly over and north of the cyclone 
N1 remains basically meridional but shifts somewhat to the 
northwest producing stronger directional shear to t¢he north- 
west of cyclone 1. In the SH, a jet maximum greater than 
60 m sec! is centered in the base of the trough 750 no 
equatorward of cyclone S, which is located under a region of 
very weak flow at mid-layer 3. A weak closed low at mid- 
layer 3 is located approximately 700 nm northwest of cyclone 
S, almost directly over the weak closed Low at mid-layer 4. 


The flow equatorward of each cyslone continues to be zonal. 


Ee COMPARISON TO OBSERVED CYCLONES 

A superficial analySis of aodel cyclogenesis might give 
the impression that these storms are realistic in every 
aspect. A closer examination reveals them to have a number 
of unusual features. In this section, the realistic and 
anomalous features will be summarized. 

Although cyclone S is more intense, the overall 
development of both model cyclones resembles that of storms 
observed in «he North Atlantic by Petterssen, et al. (1962). 
The model cyclongenesis begins in the lower troposphere in 


the Wecinity of a baroclinic zone, under nearly zonal flow 
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aloft, and with no significant vorticity advection. The 
upper-level cold trough and the low-level cyclone develop 
Simultaneously while maintaining their position relative to 
each other. ifn addition, the primary contribution to 
development 1s provided by lower-tropospheric thermal advec- 
tion (at least during the early stages of development) and 
not upper-level vorticity advection. 

During the early stages of cyclogenesis, the lower- 
tropospheric features are very realistic. Near the surface, 
the simulated cyclones resemble the classical cyclone model 
with the cold air on the west side of the cyclone and warm 
air on the east side. Cyclones N1 and S$ are nearly symmet- 
rical for the first 24 to 48 hours of this sets However, 
the subsequent development results in a very elongated 
structure. Although asymmetrical cyclones are not uncommon, 
“he extreme elongation of thes2 cyclones ioes appear to be 
atypical. 

The lower-tropospheric thermal patterns are also real- 
istic during the first half of the period studied. The warn 
advection in the vicinity of the storm during the early 
stages of development and the warm air coincident with the 


storm center during the latter stages of simulation would 
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tend to support cyclogenesis. However, the elongated 
thickness patterns and the orientation of the thermal ridges 
in each hemisphere during the later stages are unusual. 

The Nitta and Yamamoto (1974) observational study sup- 
ports the findings of Petterssen, et al. (1962). In addi- 
tion, they noted that oceanic cyclones generally predominate 
in the lower troposphere. The development of the two 
cyclones used in this study is generally limited to the 
lower troposphere although cyclone S developes weak closed 
lows at mid- and upper-troposph2ric levels. These features, 
however, are extremely weak and do not appear until 
hour 168. It is not surprising that cyclone S developes 
to the upper troposphere, however the vertical coupling 
of this development is somewhat unusual. This feature will 
be discussed in more detail later in this section. 

The above characteristics of the nodel storms which 
appear to be realistic are basically associated with the 
lower troposphere. On the other hand, the anomalous fea- 
tures are generally associated with the vertical configur- 
ation of the cyclones and the structur2 of the mid- and 
upper-tropospheric levels of the model, mid-layer 4 and 


mid-layer 3 respectively. The development of the short 
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waves in the middle and upper troposphere appears to be 
atypical. In particular, the »srientation of these waves and 
the resulting shear are quite unusual. As noted 2arlier, 
troughs in each hemisphere eventually +¢ilt back on 
themselves creating extreme directional shear. 

The height and thermal fields at mid-layers 4 and 3 are 
nearly in phase throughout the 72 hours of this study. Fur- 
thermore, a close examination of the height fields at these 
<wo levels reveals them to be nearly identical in shape and 
vertically stacked throughout the entire period sf this 
study. This is a reflection of a self-amplification process 
and is likely related to the lack of occlusion and shaliow 
structure of the mature cyclones. In addition, th2 upper 
troughs and the cyclones maintain their positions relative 
to each other throughout the period of this study. These 
two features are consistent with th2 observations of 
Petterssen, et al. (1962); however, the distance between 
cyclone S andthe weak upper lows which develop in the SH 
after 168 hours of Simulation is unusual. The closed upper 
lows are vertically aligned but they are nearly 700 nm to 
the northwest of the surface center. By the time the upper 


lows develop, cyclone S has reached its peak intensity (975 
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mb between hours 141 and 147) at. the surface and has begun 


to weaken somewhat (to 981 mb). One might expect the 
upper-level lows to be more intense and vertically aligned 
with the surface low during th mature stage of the cyclone. 

The development of the patterns in the mid-layer 4 vor- 
ticity flelds also appears to b2 anomalous. In particular, 
the location of the vorticity maxima with respect to the 
cyclones (over cyclone N1 and equatorward of cyclone S) 
during the final 48 hours of this study is unusual. The 
extremely elongated regions of relatively uniform vorticity 
which parallel the wind flow also appear to be unusual. 

In conclusion, these model storms are characterized by 
both realistic and unusual features. The realistic features 
are generally in the lower troposphers with the atypical 
features generally in the middle and upper troposphere. The 
vertical resolution of the moiel's upper troposphere could 
possibly be a factor in the development of these unusual 
features. Although all aspects of these storms are not typ- 
ical of oceanic cyclones, they are realistic s¢nough to be 
useful in determining the mass and circulation budgets of 


numerically simulated cyclones. Quasi-Lagrangian diagnostic 
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techniques are described in the next chapter. An analysis 
of the mass budgets of these storms is presented in 


Chapter V. 
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IV. QUASI-LAGRANGIAN DIAGNOSTIC TECHNIQUE 


i. QUASI-~LAGRANGIAN DIAGNOSTICS 

A primary goal of this thesis is to describe cyclone 
development quantitatively through an analysis of mass bud- 
gets. The quasi-Lagrangian diagnostics approach provides a 
means of achieving this goal. Wash (1978) demonstrated the 
validity of the quasi-Lagrangian diagnostics approach in 
budget studies of observed and numerically-predicted 
cyclenes. Basic eat, ae diagnostics techniques are 
presented in this section with the mass budget equation pre- 
sented in the following section. A description of the 
generalized quasi-Lagrangian diagnostics budget equation and 
the budget volume coordinat2 system (Wash, 1978) are 
presented in Appendix B. 

Traditionally, budget studies have amployed a stationary 
budget volume which is positioned so that the feature being 
studied moved through or remained within the budget volume 
during the period of the study (e.g. Smith (1973), Kung and 
Baker (1975), and others). Since the horizontal extent of a 


Stationary budget volume is greater than that of a cyclone 
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scale feature, it is not possible to isolate a cyclone from 
the environment in which it is moving. In addition, budget 
Studies employing stationary budget volumes provide a quali- 
tative estimate of the large scale forcing processes 
involved in cyclone development and decay. Therefore, the 
use of a stationary budget volume is not particularly suited 
to budget studies of cyclones. 

The quasi-Lagrangian technique is ideal for budget stud- 
ies cf cyclones because the budget volume moves with the 
velocity and direction of the cyclone. A translating budget 
volume affords a method of distinguishing between the 
effects of storm translation and actual cyclone development. 
Thus a translating budget volum2 is, in fact, an instrument 
for examining the fluxes across the pada ot volume boundaries 
by eliminating the distortion that results from «he movement 
of the storm. 

Budget studies which employ a storn-following budget 
volume provide a view of the interaction between the smaller 
scale cyclone and the planetary scale environment in which 
lit is embedded. Through the transport relation and other 
governing equations, the quasl-Lagrangian budget approach 


provides a method of calculating the underlying scale 
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interaction processes. Eis approach emphasizes the 
transports of mass through the boundaries and the changes in 
the circulation within the budget volume. Quasi-Lagrangian 
diagnostics incorporates transports ints and out of the 
budget volume, as well as internal ssurces and sinks of 
these basic properties. 

The guasi-Lagrangian diagnostics approach has been used 
in avariety of cyclone studies. Panola i, partially 
extracted from Wash (1978), is an updated list of severai 
recent budget studies which applied quasi-Lagrangian 


diagnostics to observed and/or simulated cyclones. 


TABLE III 
QUASI-LAGRANGIAN DIAGNOSTICS BUDGET STUDIES 
(partially extracted from Wash, 1978) 

PROPERTY RESEARCHER (S) 
Available Potential Energy Spaete 1974 
Circulation Wash 1975 
Mass Johnson and Downey 1976 
Absolute Angular Momentun Johnson and Downey 1976 
Kinetic Energy Chen and Bosart 1977 
Absolute Angular Momentun Wash 1978 
Mass Roman 1981 
Carella tion/ 

Angular Momentum Conant 1982 
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As mentioned earlier, quasi-Lagrangian transport theory 
provides a diagnostic framework which is ideally suited to 
the study of cyclone-scale circulations. In addition, the 
internally consistent "data" base of numerically simulated 
cycicnes limits computational errors to those resulting from 
Interpolation from a rectangular grid to the circular grid 
of the budget volume and to round-off errors. With the 
quasi-Lagrangian diagnostics approach, cylindrical budget 
volume is placed around the model storm and is moved with 
the storm center. Por this study, the surface pressure 
Minimum is used to indicate thea storm center. The budget 
volume extends from the earth's surface to the top of the 
atmosphere with the axis of the budget volume normal to the 
earth's surface. The budget volume is actually a portion of 
a cone with its vertex at the earth's center. However, 
Since the tropospheric height is much less than the radius 
of the earth, the budget volume radius is assumed to be 
constant in the vertical. The budget volume for this study 
consists of the seven model sigma levels as defined in 
Chapter II. The radius of the budget volume is defined in 
degrees latitude relative to the storm center. Having 


defined the location and radius of th2 budget volume, 


73 





the budget of any basic storm properties (e.g. mass, 
circulation, angular momentum, kinetic energy, and moisture) 


can te studied. 


Be EUDGET EQUATIONS 

For this study the quasi-Lagrangian jiiagnostics techni- 
que is applied to mass budgets of two numerically simulated 
cyclones. The model experiments defined vertical motion to 
be zero at the upper and lower boundaries of the atmosphere. 
Thus vertical transports across these boundaries are 
eliminated in the mass budget equations. 

Table IV contains the mass budget squation in the sigma 
coordinates of the model. In its most basic form, the gen- 
eralized budget equation consists of thr2e terms: lateral 
transport, vertical transport, and source (or sink). 

Since there are no sources or sinks of mass, the source 
(sink) term cf the generalized equation is zero. The 
remaining terms in the mass budget equation represent the 
lateral mass transport (convergence/divergence) and vertical 
transport (residual) respectively. When the limits of inte- 
gration are 9 = 0 and 09 = 1 (where og = 0), the budget resid- 
ual term must be zero. feenoc, the Lateral mass transport 


term for the various model layers are adjusted so that the 


74 






integrated residual vanishes. Lateral mass transport and 
mass tendency are computed within the mass budget and resi-~ 
dual imbalances are then computed from the time averages of 
these quantities. Vertical integration of the layer resid- 
uals produces a vertical mass transport which is then 


expressed as a vertical velocity field. 


TABLE IV 
QUASI-LAGRANGIAN MASS BUDGET EQUATION IN SIGMA COORDINATES 


(after Wash, 1978) 


(See Appensix B for budget notation 
and generalized budget equation.) 


For the mass budget f =-1 and = = O65 





Substitution of these values into the 
generalized budget 2quation results in 
the following mass budget equations 
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INTRODUCTION 

The analysis of the quasi-Lagrangian mass budget begins 
With a presentation of the thermal structure of the cyclones 
in terms of potential temperature (Section B). This is 
follewed by a discussion of the pressure (mass) tendency 
(Section C) and the teta 1 lateral mass transport 
(Section D). Vertical mass transports are discussed 
relative to the model vertical velocities (Section E). In 
an effort to determine the accuracy of the vertical velocit- 
22s computed in the mass budget, budget vertical velocities 
are compared to the model vertical velocities (Section F). 

Most figures presented in this chapter are layer/tine 
pericd cross-sections. Budget calculation values at levels 
and "data" times were averaged to produce layer and time 
period midpoint values respectively for the time ssctions. 
Exceptions are indicated wher appropriate. Inward and 
upward transports are represented by positive values in all 


these figures. 
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Be. THERMAL STRUCTURE 


Area-avertaged mid-layer and "data" time potential temp- 


eratures are presented in Figs. 28 and 29. The vertical 
structures are similar in the two csyclone2s. The patterns 
for radii six and nine are nearly identical. In the upper 


tropcsphere, the potential temperature increases slightly, 
while decreases are seen in the midile and lower troposphere 
during the last half of the study. During the period of 
most intense cyclogenesis (hours 96 to 125) only small fluc- 
tuations are seen in ‘the area-averaged potential tempera- 
ture. The cooling, which occurs during a quasi-steady state 
pericd for each cyclone, iS more pronounced in cyclone N1. 
This cooling within each budget volume is a consequence of 
the development and poleward movement of the cyclones. In 


addition, although the troposphere of syclone N1 cools more 


‘ham cyclone S, its average potential temperature is orly 
Slightly lower at the final "data" time, Since cyclone N1 
was initially warmer in the low troposphere. These 


potential temperature cross-sections show that cyclone N1 
initially had weaker vertical stability than cyclone S. 
However, by the end of the study, the two cyclones have 


Similar vertical stability characteristics. 
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The cooling in cyclone N1 gives the appearance of a cold 
core structure. This however is not the case since the cen- 
ter of each cyclone is actually located in the warmest air 
Within each budget volume (see Fig. 11). The cooling is 
likely a result of cold advection into the equatorward half 
of each budget volume. The low-tropospheric thermal pattern 
can be seen in the 1000 mb to mid-layer 4% thickness pattern 
fer hour 168 (see Fig. 11). The thickn2ss pattern shows a 
warm tongue extending from the east through the center of 
cyclone N1, while a similar warm tongue extends east to west 
through the radius nine budget volume of cyclone S. The 
thickness pattern also reveals that the cold air occupies 
more of the budget volume of cyclone N1 than cyclone S, thus 
resulting in the slightly lower average temeperatures for 
cyclicne Ni. 

Although the model atmosphere has nd specified tropo- 
pause, Figs. 28 and 29 clearly show the existence of a 
strong tropopause and a stratosphere during the 72 hours of 
this simulation study. Level 2 (approximately 200 mb) 
effectively functions similar to a ‘tropopause layer with 
slight warming and strong vertical stability above 200 mb. 


The potential temperature evolutions ar2 characterized by 
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Slight fluctuations below level 2. These fluctuations 
provide a hint of a cyclic evolution pattern during the 
development of these storms, although the period of the 


cycle is not evident. 


C. MASS TENDENCY 

There are no sources or sinks of mass within the atmos- 
phere, however it is possible for a limit2i volume of atmos- 
phere to gain or lose mass. The time rat2 of change of mass 
within a storm-following budget volume provides a quantita- 
tive measure of cyclogenesis. Within th2 mass budget, the 
total mass within each specified budget volume is computed 
for each "data" time. From these values, the mass tendency 
between "data" times is determined, and subsequently 
converted to a pressure tendency. 

The mass tendency at 3-hour intervals for both cyclones 
iS presénted in Fig. 30. A cyclic pattern is clearly 
indicated. These time profil2s also indicate that the 
development of the model cyclones occurs in two phases--an 
initial explosive cyclogenesis period (hours 96 to 126) and 
a quasi-steady state period (hours 126 to 168). A 12-hour 
oscillation is clearly evident throughout these two phases. 


During the explosive cyclogenesis phase, the cycle in the 
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radius six budget consists of a maximum 3-hour pressure fall 
followed by a 9-hour period of relatively smaller surface 
pressure change. The cyclic pattern during the quasi-steady 
state phase is characterized by alternating periods of 
deepening and filling. The 12-hour cycle can be seen in the 
pressure tendency patterns for both th2 inner and outer 
radii. The pressure tendency trace for each budget volume 
shows that thecyclone's area-averaged surface pressure drops 
in discrete intervals like a stap function rather than show- 
ing continuous deepening. These oscillations could be 
removed by using an highpass filter. Filtering would also 
serve to emphasize the deepening during the first 24 hours 
eo “ne Tiliing Which occurs later. 

The pressure tendency oscillations ar? “on as pronounced 
in the graphs of the central surface pressures of these 
eyecleones w(Fig. 31). This is especially true during the 
explcsive cyclogenesis phase, when a generally steady 
decrease is seen in the surface pressur2. However, the rate 
of change nearly corresponds to the pressure tendency fluc- 
E@@tions in Fig. 30. During the. second phase of cyclone 
development, oscillations are seen in the central surface 
pressure, however they are not as large as those the 


pressure tendency. 
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The cause of these mass fluctuations is not known. One 
possibility is that the cycle is the reflaction of an exter- 
nal gravity wave propagating within the model. An external 
wave would have a Similar effect on the surface pressures of 
neighboring anticyclones. For this reason a study of the 
central surface pressures of mcedel-generated anticyclones 
was completed. 

This preliminary study of five highs (three in the NH 
and two in the SH) for periods of 48 to 72 hours (the length 
of the period was dependent on whether the anticyclone was 
an analyzed feature at a particular "data" time) strongly 
indicates a cyclic pattern of surface pressure fluctuations 
(approximately two millibars) in the model anticyclones as 
well. The pressure oscillations of the cyclones and anticy- 
clones suggest a propagating external wave within the model. 

The pressure tendency traces show larger fluctuations 
for radius nine than radius six. This indicates that the 
pressure falls/rises for the area between radii six and nine 
are larger than those for the area within the radius six 
budget circle. An examination of other budget radii showed 
that the magnitude of the pressure teniency fluctuations 


znereases for larger Badia. This also suggests a 
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large-scale (perhaps hemispheric or global) Sloshing. 
Because the scale of this feature is nuch larger, and it 
does not interact with the synoptic-scale features, the 
Slosh is independent of cyclone development. A’ more 
detailed investigation of the model surface pressure field 
1s required to isolate this phenomena and determine if it 


adversely affects the numerical simulations. 


D. LATERAL MASS TRANSPORT 

Cne goal of this thesis is to isolate the lateral and 
vertical branches of the mass circulation of the model 
cyclones. The lateral mass transport which is the mass flux 
through a specified budget volume boundary. The lateral 
mass transport time sections (Figs. 32 and 33) graphically 
describe a classical cyclone structure of low-level converg- 
ence (surface *to approximately level 4: surface to 600 mb) 
and up per-level divergence (approximately level 4 to 
mid-layer 2: 600 to 150 mb). A Similar pattern is noted at 
radius nine, however the level of nondivergence (LND) is 
slightly higher (approximately nid-layer 4: 500 mb). The 
lateral mass transport time sections indicate that the mass 
circulation does not extend into the top two model layers. 
This further substantiates the existance of tropopause and 


stratosphere behavior within the model. 
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The maximum inflow at the inner and outer radii of each 
cyclone occurs in the lowest model layer (centered at 
apprcximately 900 mb), while the maximum outflow is centered 
near mid-layer 3 (approximately 300 mb) at the tropopause. 
i @wddition, the mass circulation at radius nine of each 
cyclone is more intense during the quasi-steady state phase 
than during the rapid cyclogenasis phase. This is a further 
sndication that the the second phas2? of development is actu- 
ally an outward growth period more than a deepening period 
for these cyclones. 

A cyclic pattern is evident at the Levels of maximun 
inflow and outflow at both the inner and outer radii of each 
cyclone. Looking first at cyclone N1 (Fig. 32), the maxinun 
outflow at radius six oocccurs at hour 117 with secondary 
Maxima centered at hours 1480, 150, and 168 while nearly 
equal inflow maxima are seen 2t approximately these times. 
The maxima at hour 117 occur iuring the period of intense 
surface development (as indicated by changes in the central 
surface pressure, Fig. 31), while the maxima at later times 
eccur ijiuring the quasi-steady state phase of cyclone Nl. 
The lateral ‘mass transports at radius nine of cyclone N1 


are also characterized by oscillations in ‘ths basic 
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ant lom-outftlew pattern. Tye Gyelic® @ariation in the 
lateral mass transport are coupled with the osciliations 


in the pressure field. 
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A dramatic increase at 69 and a decrease at 99 in the 
height of the level of LND is also seen during the quasi- 
steady state phase of cyclone N1. The LND height change for 
radius nine (approximately hour 135) is related to the 
intensification of the mass circulation during the second 
stage of development and the growth of the cyclone to the 99 
budge~ volume boundary. 

The lateral transport time sections for cyclone S also 


exhibit the classical mass circulation pattern. However, 
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*he cyclone S lateral transport pattern is noticeably dif- 
ferent from that of cyclone Ni. Specifically, the cycione S 
lateral mass transports at radius six ar2 more intense «han 
those of cyclone N1. This is a reflection of the greater 
intensity of the spring hemisphere cyclone. In addition, 
the height of ‘the LND at radius six of cyclone S gradually 
lowers from about 500 mb to near 700 ob, and is generally 
lower than its NH counterpart (which remains near 600 mb 
throughout this study) after hour 132. 

The radius nine time sections for these *wo cyclones 
also display similarities and differences. As in 
cyclene N1, there is a Simultaneous decrease in the height 
of the END and the intensification of the lateral mass 
transports. However, this occurs earlier in cyclone S (hour 
w20 for cyclone S and hour 135 for cyclone Nt), and the 
decrease in «he LND height is ore dramatic (from 500 mb to 
700 mb for cyclone S and about 475 ab to 600 mb for 
cyclone N1). The sharp decreas? in the height of the radius 
nine LND in cyclone S occurs simultaneously with the outward 
growth of the cyclone to the outer budget volume boundary. 
AS with radius Six, the radius nine lateral transports are 


more intense in cyclone S. 
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Of special note is the abrupt weakening of the cyclone s 
mass circulation during the final -21 hours of this study. 
This weakening is most noticeable and occurs earlier at 
radius six. The decrease in the mass transport intensity 
can te seen in both the inflow and outflow branches, however 
it is more pronounced in the outflow branch. This waakening 
of the lateral mass transports is indicative of the decay of 
cyclone S and is reflected in a correspsnding increase in 
the central surface pressure of cyclone S. In contrast, the 
lateral mass transports at the cyclone N1 budget radii con- 
tinues to be strong at the last "data" time when cyclone N1 


reaches its lowest central surface pressure (984 mb). 


Bie VERTICAL MASS TRANSPORT 

The vertical branch is an integral part of a cyclone's 
agsoe circulmtion. Since the vertical mass transports are 
the area-averaged vertical velocity structures within the 
cyclene (in both isobaric and sigma coordinates), time sec- 
tions of the model vertical velocities (Figs. 34 and 35) are 
used for the following discussion of the vertical rans- 
ports. This discussion vooints out Similarities and differ- 
ences in the patterns of thes= two storms. In general, 


upward motion is seen below level 2 (approximately 200 mb) 
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at the inner and outer radii. Alternating periods of weak 
upward and downward motion are seen above level 2 in the 
model stratosphere. The upward motion is quite intense at 
radius six of each cyclone. The level of maxinum upward 
vertical motion coincides with the LND indicated in the lat- 
eral mass transport time séctions (Figs. 32 and 33). 
periods of maximum upward vertical motion coincied (as they 
must) with the periods of most intense lateral mass trans- 
port. transport. The upward vertical motion at radius nine 
is weaker than at radius six. However, the intense upward 
vertical motion occurs only within the troposphere. 

Considering first the six degree budget volume, two dis- 
tinct periods of maximum upward vertical motion are evident 
in each cyclone. For cyclone N1, the first period (hour 107 
to 123) occurs during the explosive cyclogenesis phase, 
while the second period (after hour 147) occurs during the 
quasi-steady state phase. The subsequent radius six maximum 
occurs simultaneously with a similar maximum at radius nine, 
and results from the outward growth of the cyclone to the 
outer budget volume boundary. 

Cyclone S$ also experiences two periods of intense upward 


vertical motion (centered at hours 117 and 138). As in 
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cyclone N1, the first maxinum occurs during the explosive 
cyclogenesis phase while the second maxiaum and its radius 
nine counterpart occur during the period of outward lateral 
g@row#th. Near ‘the beginning of the integration, qa cneEa 
period sf intense upward vertical motion is evident for each 
budget radius. dowever, its significance is unclear due to 
the lack of the preceeding nodel data. In addition, the 
decay of the cyclone is indicated by the rapid decrease in 
upward vertical motion after hour 159. 

An interesting feature of the cyclone S radius nine +ine 
section is the downward vertical motion indicated in the 
lower t~wo layers between hours 101 and 111. This requires a 
downward motion within the 6°9-9° area which more than off- 
sets the upward motion within radius six. This indicates 
strong subsidence around the cyclone compensating for the 
intense upward vertical motion associated with the explosive 
cyclogenesis phase when cyclon2 S$ is still within the 69 
budget volume. This is a possible reflection of the cyclic 
evolution of the model cyclones. 

Time traces of area-averaged vertizsal velocities for 
level 4 of cyclones N1 and S ar2 presented in Fig. 36. The 


+wo major periods of strong upward vertical motion which 
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correspond to the two-phase cyclone development are clearly 
depicted. The outward growth phase of each cyclone is viv- 
idly displayed in these traces. A dramatic increase in the 
radius nine area-averaged vertical velocities is seen as the 
lateral growth of these cyclones exceeds the 9° budget 
volume boundary. 

The radius six trace for cyclone N is characterized by 
two distinct maxima at hours 115 and 168. The radius nine 
trace essentially has a singl® maximum which occurs during 
the second phase of cyclone development. The extreme max- 
imum at hour 115 iS approximately in the niddle of the rapid 
cyclogenesis phase when the central surface pressure 
decreases 17 mb between hours 138 and 123 (see Fig. 31). 

A surprising feature of the vertical velocities occurs 
between hours 141 and 150 when the magnitude within radius 
nine actually exceeds the corresponding value for radius 
SUX. This implies that the vertical velocities within the 
69-99 region exceed those within the radius six budget 
Grr che. 

Although the cyclic pattern seen in cyclone S is similar 
to that of cyclone N, there ar2 significant differences. As 


im cyclone N, two maxima are seen at radius six 
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(hours 118 and 136), however the time interval between these 
two cccurrances is much shontar than in cyclone N. The 
first maximum occurs during the first phase of cyclone 
development when the central surface pressure decreases 
22 mb between hours 99 and 126. The second maximum, which 
occurs during the quasi-steady state phase of cyclone S, 
coincides with the single maximum at radius nine. The maxi- 
mum at radius nine of cyclone S is larger than its cyclone 
N counterpart. The decay of cyclone S is clearly indicated 
by «he sharp decrease in the level 4 vertical velocities at 


both six and nine degrees. 


e. COMPARISON OF MODEL AND BUDGET VERTICAL VELOCITIES 

One objective of this thesis is a comparison of the 
model vertical velocities with those computed in the mass 
budget. The vertical divergence of the vaetieall mass trans- 
port (budget residual) and the mass tendency are computed 
within the mass budget (see Tables IV and V). A vertical 
integration of these values (from an upper boundary value of 
= © a = 0) describes the vertical mass transport, which 
is then converted to an area-averaged vertical velocity. 
Thus the budget vertical velocities are essentially kine- 


matic omega values for the specified budget volume. 
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Figuze 36. Model vertical veloc for level 4. 
(- “4 bar sec!) as apa bp. aS in Fig. 30. 
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This approach has been used in atmospheric studies to 
infer vertical velocities and vertical transports of mass 
and other properties. In this study, vertical velocities 
ac2 available as part of the model "data" base. This 
affords the unique opportunity to test the accuracy of the 
budget-derived vertical velocities. 

The basic patterns in the aodel vertical aaa time 
sections (Pigs. 34 and 35) are contained in the budget ver- 
tical velocity time sections (Pigs. 37 ani 38). The budget 
vertical velocities for radii six and nine of each cyclone 
are most intense below the tropopause (approximately nid- 
bager 3). Very weak vertical motion is seen in the model 
stratosphere. At radius six of each cyclone, two maxima are 
Seg in the budget vertical velocities time sections, with 
the earlier maximum occurring during the first phase of 
cyclone development. AS in the model vertical motion pat- 
tern, the second maximum in the budget vertical velocity 
time section occurs in conjunction with a radius nine maxi- 
mum. This again indicates the lateral growth of the 
cyclones beyond «he 9° budget volume boundary. For each 
budget radius of cyclone S there is a significant decrease 


in the budget vertical velocities after hour 159 when 
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G¥Yctome S begins to fill. This feature is also seen in the 
model vertical velocities. 

Although the patterns in the model and budget ucees 
velocity time sections are similar, significant differences 
were noted in the magnitudes of the budget vertical veloci- 
ties. Timesections of the differences between these two 
values (budget value minus mod21l value) are presented in 
Pigs. 39 and 40. With the exception of the stratosphere, 
the budget vertical velocities were larger than those of the 
model. It appears that the largest departures occur where 
the vertical velocities (both model and budget) are largest. 

Since the model has an internally consistent "data" base 
there are two possible sources for the errors in the budget- 
derived vertical velocities. Errors can result from the 
interpolation of the winds from the rectangular grii of the 
model to the circular grid of the mass budget. Other errors 
arise from the vertical distribution of the imbalances 
between the integrated divergence and the mass tendency, 
which must balance since there are no sources or sinks of 
mass. The comparison of the budget and model vertical 
velocities indicates that the lateral mass transports are 


“oo large after the interpolation and balancing proceedures 
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are completed within the mass budget. Given the consistent 
"data" of the model, these errors are generally larger than 
one might anticipate. A detailed examination of the inter- 
polation scheme and the vertical distribution of the budget 
residuals should be conducted in an effort to isolate the 


source of these errors. 
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VI. CONCLUSIONS 


In this thesis two numerically simulated oceanic extra- 
tropical cyclones were studied. Initially a synoptic study 
of these cyclones, including a comparison with observed 
cyclones, was completed. These cyclones were found to be 
characterized by both realistic and unusual synoptic fea- 
tures. In a general sense, the realistic features were 
associated with the lower troposphere while the unusual fea- 
tures were associated with the middle and upper troposphere 
and the vertical coupling of the model layers. 

These cyclones also were studied through the application 
of quasi-Lagrangian diagnostics. This diagnostic study of 
the mass circulation of model cyclones showed a concentra- 
*ion of inward lateral transport (convergence) in the lowest 
model layer and outward lateral transport (divergence) ina 
layer centered near 300 mb. Although the model has no spec- 
ified tropopause, the basic patterns o£ the lateral and 
vertical mass transports and the thermal structure revealed 
«he existence of a tropopause near mid-iayer 3 of the model. 

Traces of the mass tendency aera 69 and 9° budget 


volumes revealed a 12-hour cyclic pattern imbedded in a 
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two-phase cyclogenesis pattern. The first phase of cyclone 
development was a period of explosive daepening while the 
second phase was a quasi-steady state period. The 12-hour 
cycle was most easily seen in the pressure tendency trace. 
During the first phase of cyclone development, the cycle 
consisted of a maximum 3-hour pressure fall followed bya 
9-hour period of little surface pressure change. During the 
second phase, the cycle was characterized by alternating 
pericds of pressure falls and rises. It was also noted that 
the fluctuations in the pressur2 tendency incr2ased for suc- 
cessively larger budget radii. The two-phas2 cyclone 
development was seen in the traces of the central surface 
pressures. However, the 12-hour cyclic pattern was not eas- 
ily detected in these traces. This 12-hour cycle produces a 
step function decrease in the area-averaged surface pressure 
of these cyclones rather than a steady deepening. This step 
function evolution is not clearly indicated by the central 
surface pressure traces (see Fig. 31). A possible cause of 
this cycle is an external wave propagating within the model. 
However, further studies are naaded to isolate the source of 


these oscillations. 
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The lateral mass transport and the vertical mass trans- 
port (converted ‘0 a vertical velocity) clearly described 
the mass circulation of these cyclones. Time sections of 
these terms also indicated two-phase cyclone develpment. 
TwO maxima were seen at radius six of each cyclone--one 
during each phase of development. The ne maximum seen at 
radius nine occurred in conjunction with the second maximum 
at radius six and indicated the outward growth of the 
cyclones beyond the outer budget volume boundary. Maxima in 
*he lateral mass ‘transport and vertical velocity time sec- 
tions generally occurred simuita neously. In addition, the 
onset of the decay of cyclone S was indicated by the rapid 
decrease in the intensity of the lateral mass tendency and 
vertical motion near the end of the study. . 

Area-averaged vertical velocities were computed within 
the mass budget and were compared to area-averaged model 
vertical velocities. Although tha time sections of these 
terms have similar basic patterns, significant differences 
were found in their magnitudes. It appeared that the depar- 
+ures were largest where the vertical velocity (both model 
and budget) were the largest. Th= two seget nie sources of 


thes® errors are interpolation of the wind from ‘the model 
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grid to the budget volume grid, and imbalances in the 
vertical distribution of the mass budget residuals. Further 
studies are required to isolate the relative magnitudes of 


these two errors. 
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VII. RECOMMENDATIONS FOR FURTHER STUDIES 


The results of the synoptic study and the mass budget analy- 
Sis suggest that further diagnostic studies of the basic 
properties of these cyclones are needed to fully understand 
+heir dynamic structure. In addition, further studies are 
needed to determine the source of the errors in the budget 
derived vertical velocities. With this in mind, the follow- 
ing recommendations are made: 


e Examine the budget technique for computing area~ 
averaged vertical velocities and the interpolation 
Scheme employed within the mass budget to isolate the 
4 ae 9f the errors in the budgéd derived vertical 
velocities. 


e Perform budget studies on the model cyclones. duri 
eo and later time p2riods than covered in th 
study. 


ng 
Ls 


e Perform budget studies on model generated anticyclones 
to determine if ‘they are characterized by the ?2-hour 
eyelve pettern of the model cyclones. 


e Study the CYGwic.. oscillations by differencing 
successive pressure fields. 


e Perform budget studies on See cones 
generated in the Sandgathe (19 
experiments. 


yclones 
a: 


and antic 
8 1) adiabatic 


e Parform budget studies of other ovroperties on these 
model cyclones and on model anticyclones. 
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RPPENDIX A 


COMMONLY OSED ACRONYMS AND ABBREVIATIONS 


NOS APS 


UCLA 


NNW 


ENE 


SSE 


NNE 


SSW 


LND 


AVA 


CVA 


PBL 


sec~? 


mb 


me 


OK 


Navy Operational Atmospheric Prediction Model 


University of California, Los Angeles 


Naval Environmental Prediction Research 
Paciiity 


east-northeast 


east-southeast 


north-northeast 


south-southwest 


level of nondiverg ence 


anticyclonic vorticity advection 


cyclonic vorticity advection 


planetary boundary layer 


per second 


eidilibar 


degree Celcius 


degree Kelvin 








NA 


SH 


na 


northern eee IS ST 
southern hemisphere 
nautical miles 

mass 

north 

south 


meters 





APPENDIX B 


QUASI-LAGRANGIAN BUDGET FRAMEWORK AND 
GENERALIZED BUDGET EQUATIONS 


The purpose of this appendix is to present the formula- 
tion of the storm budget volume coordinate system and the 
generalized budget equations which can be applied to any 
desired atmospheric property. The information contained 
herein is extracted from Wash (1978) as adopted from Johnson 
and Downey (1975 a). 

The storm budget volume in spherical coordinates is 
@epaicted in Fig. #1. The following notation is relevant to 


this figure: 


of Azimuthal coordinawe 
B Angular radial coordinate 
g Position vector from 2arth's 


center 


R Poeutaon Secetor frou storm center 
= to any point in budget volume 
m,1l,k Orthognal unit vectors along 


meridional, azimuthal and ; 
vertical directions of spherical 
coordinates for storm 


Wind velocity relative to earth 


Vea 


lates 


— 





{x 


Horizontal velocity of reference 
axis 


Horizontal Hee of any point 
on the volume boundary. 


Table VI contains a listing of the generalized budget 


eguations. The following is an explanation of the notation 


used in these equations: 


-W 
C07 Hi, 


(dg/dt) 


csinp 


'~ 


th 


The vertical coordinate systen 


Volume 


Top boundary 


Bottom boundary 


Horizontal wind component relative 
to budget volume normal to lateral 
boundary 


Vertical velocity divergence into 
a layer within the volume 


Radius of budget volume 


Surface pressure 


Specific budget property (per unit amass) 


Vertical velocity (dao /dt) 


Acceleration dus to gravity 











horrzontal 
surfaces 


Eartn’s 
surface 


5 


Figure 41. eae emerangsan Storm budget volume coordinate 
System. 


115 





. TABLE V 
GENERALIZED QUASI-LAGRANGIAN BUDGET EQUATIONS IN SIGMA 
COORDINATES 
(after Wash, 1978) 


The budget volume integral is defined as 


27 

Or Ba — | 

F = | | ‘+ fr~ sinfdod8(-do). 
ao 


The most general form of the budget equation is 


dF 


= * LT (F) + #£4VI(F) + ~~ S(F) 


where the lateral transport term is 
LT(F) J [2 — = £(U-) gr sin Bda(-do) |8, 


the vertical transport tera is 


Or Be oT 
bd P j 
VT (F) 7 * 9” (£0) sin8dad8 (-do), 
O. 0 0 g a0 


and the source (sink) tera is 
O 


Ba y 
5 (F) =f J Ll? ee a’ sin8dad8(-do) . 


Oo; 


The vertical transport term is zero 
when the limits of integration are 
;* 1 and of = 0,since there can be 
transport through the top and 
bottor poundaries of the budget volume. 
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